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p ■ Abstract 

We show that a simple scaling model of very forward particle production, 

consistent with accelerator and air shower data, can describe all features of 

the very high-energy interactions recorded with emulsion chambers. This 

is somewhat surprising after numerous claims that the same data implied 

large scaling violations or new dynamics. Interestingly, we cannot describe 

some of the Centauro events, suggesting that these events are anomalous 

independently of their well-advertised unusual features such as the absence of 

neutral secondaries. 
PACS number (s): 13.85 
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I. INTRODUCTION 



In this paper we show that the simplest assumptions on particle production in the very 
forward direction, Feynman scaling and constant inelasticity, can provide a coherent descrip- 
tion of air shower and emulsion data. All features of our model are approximately shared 
with QCD-inspired models |I[ of high-energy particle interactions. The phase space popu- 
lated by the secondaries in very- high-energy cosmic-ray interactions is, unfortunately, not 
covered by collider experiments which only view centrally produced particles. We therefore 
start the discussion with the simplest ansatz for the dynamics of particle production at large 
rapidities. In the end, this guess will turn out to be adequate in describing a wide variety 
of data. 

In order to describe the underlying high-energy particle production dynamics common to 
accelerator, air shower and emulsion chamber data, we start with the simplest phenomeno- 
logical model for the rapidity density distribution of charged pion production which we 
parametrize as 

dN dN (l-x) n 

— — = x—— = a , (1) 

dy dx x m 

where y is the rapidity of the produced secondaries. The Feynman variable x is given by 
the ratio of the secondary particle energy E to the incident energy E . 

The (1 — x) power law in (|]) describes the decrease in the number of particles produced 
with large x values, the 1/x term dictates the shape of the function at low x. Simple parton 
counting rules [0 imply that n ~3 and m = — 1/2. Eq. ([J) reflects the approximate 
Feynman scaling behavior in the fragmentation region (large x, y) expected from QCD [|TJ. 
Our basic assumptions imply that the inelasticity K, i.e. the fraction of the collision energy 
going into the production of secondaries, is independent of energy. QCD-inspired models 
predict that K = 0.6 — 0.65 [|3|-p at the energy region of interest. We have to account for 
nuclear target effects which can simply be done by allowing increased m-values to describe 
the enhanced production of low- x particles in nuclear collisions. Also, nuclear effects increase 
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the value of the proton-air inelasticity || to approximately 0.8. This simple picture is consis- 
tent with the little information one can extract from accelerator results on the production of 
particles in the very forward direction |7j . Especially in the high-energy collider experiments 
the phase space relevant to cosmic-ray data is usually obscured by the beam-pipe ||. 
With this background we may investigate air shower and emulsion chamber results. 



II. FORWARD PARTICLE PRODUCTION AND AIR SHOWERS 

A QCD-based Monte Carlo simulation of ultra-high-energy cosmic-ray hadron interac- 
tions, which incorporates large p t jet production, has been used by L.K. Ding et al. |9] 
(hereafter referred to as DKTY) to perform a detailed description of air shower observa- 
tions. The rapidity distribution of secondary particles has been derived in this analysis |L0 



e.g. for interaction energies of 10 16 eV (^/i ~ 4.3 TeV) as illustrated in Fig. |IV|. Also 
shown in the figure is an adequate description of the resulting distribution using Eq. ([[]) 
with a = 0.12, n = 2.6 and m — 1, parameters quite consistent with our expectations. The 
analysis of reference also shows that particle production in the forward direction exhibits 
approximate Feynman scaling, consistent with our simple assumptions. 

The main point of this paper is that this model, in its simplest form, describes all the 
details of emulsion chamber experiments. This was somewhat surprising to us given the 
many claims of new dynamics and non-scaling behavior in the literature based on the same 
information; see e.g. Ref. pT| . Interestingly, we cannot describe some of the Centauro 



events. The rapidity distribution of the secondaries in these events pegs them as being 
anomalous, independently of their well-advertised unusual features such as the absence of 
neutral secondaries. 



III. EMULSION CHAMBER DATA 



Before proceeding with the study of emulsion chamber data, we first comment on a 
crucial aspect of the calculation: the mean inelasticity which is calculated through energy 
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conservation (using charged secondaries only) 

§Wft = (2) 

The unphysical singularity of Eq. ([[]) on x = is removed by the threshold effects in the 
measurement. The lowest energy particles in the cascade will not make it to observation 
level. We implement this by making the x-distribution constant below a certain value Xq, 
so that the mean inelasticity is 



(K) = \a 



1 (f ~x) n (1 + zgT 
ax + _ x 



(3) 



The existence of a lower limit on x, resulting from inevitable energy-loss in the air cascade 



before it reaches the detector, is well known |12| . Emulsion chambers observe the debris of 
air showers initiated by cosmic rays one or more nucleonic mean-free-paths (mfp) above 
the detector. From the study of Ref. [|I2[] one concludes that, because of the strong energy 
dissipation in the cascade, Xq must be around 1 x 10~ 3 and 3 x 10~ 2 for mountain altitude 
measurements. We can obtain the same result from a back-of-the-envelope calculation. 



Take as an example the Ursa Maior event observed by the Brazil- Japan Collaboration [13 



Triangulation of the secondary charged tracks predicts that the first interaction occurred 2 or 
3 km above the chamber, which represents approximately 2 mfp. At the first interaction the 
nucleon with energy Eq releases on average around \Eq to particle production and therefore 
\Eq to each pionic component (it , 7r ± ) . After 2 more interactions, a charged pion would 
have on average (Eq/Q) 3 so that Xq ~ (1/6) 3 ~ 5 x 10~ 3 . This value is consistent with the 
one obtained in Ref. |12[ and, if used in Eq. @, leads to (K) = 0.81, matching the expected 
value; see Section |I[ 

We will investigate experimental data obtained in large emulsion chambers, using rigorous 
analytical solutions of the diffusion equations for the hadronic cascade induced by one single 
nucleon in the atmosphere. The one-dimensional solution F H (E, E ,t,t ) enables us to 



calculate the integral energy spectrum of hadrons [|T4| 



I H (> E, E ,t f , t ) = / dE'F H (E', E ,t f , t ) , (4) 

J E 



which gives the number of hadrons with energy above E, at the detection level tf for a cascade 
initiated by a nucleon with energy E at atmospheric depth t - The three-dimensional 
solution Fh(E, E , r, t, t ) is used to calculate the energy-weighted lateral spread |15],[16] 

N H (> Er, E , t f , to) = J J dE'dr' H(E'r' - Er) 

xF H (E',E ,r',t f ,t ) , (5) 

which evaluates the number of hadrons with energy times radial length greater than Er. 
The function H(E'r' — Er) is the Heaviside function. 

The only free parameters are the primary particle energy Eq and the atmospheric depth 
to of the first interaction. In our calculations we define the visible energy E^ = k 1 E 
(with k 1 ~ (k 7 ) = 0.25) and the atmospheric depths t in units of nucleonic mfp, with 
\n = 80 g/cm 2 . We use for the transverse momentum (p t ) = 0.4 GeV/c. 



In Fig. IV we compare our results for the integral spectrum of the event Ursa Maior [13 



detected at Mt. Chacaltaya (540 g/cm 2 or t f = 6.75). We obtain E^ ] = 2930 TeV, which 
corresponds in the center of mass system (c.m.s.) to yfs = 4.7 TeV. The determined depth 
of the first interaction (to = 4.95) is consistent with the expectation that the event started 
about 2 mfp above the chamber (see Table |). There is no data available on the energy- 
weighted lateral spread for this event. 

In Fig. fV| we show (a) the integral spectrum and (b) the energy-weighted lateral spread 



for the superfamily P3'-C5-505, detected by the Pamir Collaboration [|Tj]] at the detection 



level tf = 7.45 (or 596 g/cm 2 ). In the same figure we present the results for the event 



CentauroVII, detected by the Brazil- Japan Collaboration [|18j. For clarity, we have shifted 
in the figure the data of CentauroVII by a factor 10. The characteristics of both events are 
again summarized in Table |. 



The investigation of Centaurol, also detected by the Brazil- Japan Collaboration ||T9| , |20 
is presented in Fig. [IV] (see also Table |). We observe that neither the integral spectrum, 
Fig. [TV|(a), nor the energy- weighted lateral spread, Fig. [P7|(b), are reproduced by the an- 
alytical calculation. We verified that the origin of this discrepancy cannot be associated 



with an increase of the transverse momentum of the secondaries, for which there is some 
experimental evidence. 



IV. CONCLUSIONS 

Guided by Q CD- inspired models, we proposed a simple phenomeno logical model for the 
production of secondaries, based on approximate scaling in the fragmentation region. The 
model parametrization is consistent with the behaviour of air showers at ultra-high energies. 
We subsequently applied the model to the description of experimental data obtained in 
emulsion chambers at different mountain altitudes. Both integral spectrum and lateral 
spread of cosmic-ray superfamilies are successfully reproduced. 

We noticed that CentauroVII is well described by the model. From our point of view, 
it is not really different from any other event resulting from normal strong interactions. 
Actually, it has been pointed out previously (15| that this event is not a "hard" Centauro, in 



the sense that it possesses a comparable number of hadrons and gamma-rays arriving at the 
top of the upper detection chamber. On the other hand, the investigation demonstrates that 
an event like Centauro I cannot be described with the same simple hypotheses, suggesting 
that this event is clearly anomalous, even before one considers their charge-to-neutral ratio. 
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TABLES 



TABLE I. Determined parameters of emulsion chamber events analyzed. 

Event (TeV) t (X N units) */-*o (Ajv units) 

UrsaMaior 2930 4.95 1.80 

P3'-C5-505 3220 1.45 6.00 

CentauroVII 7380 1.05 5.70 

Centaurol 1890 5.75 1.00 
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FIGURES 

FIG. 1. Rapidity distribution of forward secondaries as derived from air showers with energy 
near 10 16 eV, calculated by DKTY (*) and Eq. (|l]), with a = 0.12, n = 2.6 and m = 1 (solid line). 

FIG. 2. Integral energy spectrum for the Ursa Maior event (0), detected at Mt. Chacaltaya [13], 
compared to the analytical calculation (solid line) using the x-distribution of Eq. (|l|). The fitted 
parameters and to are given in Table |[ 

FIG. 3. Description of the superfamily P3'-C5-505 (0)> detected by the Pamir Collabora- 



tion [17], and of the event CentauroVII (A), detected by the Brazil-Japan Collaboration [^q| . For 
illustrative purposes, the data of CentauroVII have been shifted by a factor 10. Solid lines are 
the analytical calculations for: (a) Integral energy spectrum, (b) energy-weighted lateral spread; 
parameters in Table Q. 

FIG. 4. Investigation of the event Centaurol (<C>), detected at Mt. Chacaltaya |19,20fl: (a) Inte- 
gral energy spectrum, (b) energy-weighted lateral spread. Solid lines are the analytical calculations; 
parameters in Table |. 
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